Resonant polaron-assisted tunneling of strongly interacting electrons through 
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The problem of resonant transport of strongly interacting electrons through a one-dimensional 
single-level vibrating quantum dot is being considered. In this paper, we generalize the Komnik and 
Gogolin model^ for the single-electron transistor with g = 1/2- Luttinger liquid leads to the case 
of strong electron- vibron interaction in a quantum dot. The effective transmission coefficient and 
differential conductance of the system has been derived for the general case of asymmetric tunnel 
barriers. The main result obtained is that, in the zero-temperature limit, the resonant polaron- 
assisted tunneling with perfect transmission is possible. This resonant tunneling is of the novel 
(Andreev-like) type due to a special electron-electron interaction in the leads. As a result, a strong 
domination of resonant polaron-assisted electron transport at low temperatures has been found. 
Additional narrowing due to electron-electron interaction in the leads, is roughly the same for all 
polaron-assisted resonances. 

PACS numbers: 73.63.Kv, 72.10.Pm, 73.23.-b 
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I. INTRODUCTION 

The problem of electron transport through various 
molecular structures and, particularly, through molecular 
single-electron transistors (SETs)has become a hot topic 
in the modern mesoscopic physical—. A single-electron 
transistor in question is modelled as a single-level quan- 
tum dot (QD), weakly coupled to two one-dimensional 
leads (quantum wires or carbon nanotubes)^^ ' 11 ' 13 ' 14 . 
One can control the tunnel current through a de- 
scribed system by means of two independent parame- 
ters: (i) "driving" voltage (bias) V between the leads 
and, (ii) "gate" voltage V g , which is able to move the 
fcrmionic level of the QD. Besides that, the average 
current depends on the properties of tunnel barriers, 
on the electron- vibron coupling in the QD, and on the 
strength of electron-electron interaction in the system. 
For the most general case of arbitrary electron-electron 
and electron- vibron interactions, it is impossible to solve 
the transport problem exactly. Nevertheless, some limit- 
ing cases are solvable. 

It is widely knowr£~— ' 1 1 ' 16 ' 17 , that the vibrations of QD 
qualitatively change the character of low-temperature 
electron transport through a SET even with Fermi- leads. 
Two main effects characterize the low-temperature dif- 
ferential conductance of the system with vibrating QD 
in the case of Fermi leads: (i) vibron-assisted tunnel- 
ing, where the additional (satellite) resonant peaks of 
differential conductance emerge^—, and (ii) the pola- 
ronic narrowing (Frank-Condon blockade) of the widths 
of all resonances 7,16 (for strong enough electron- vibron 
coupling). 

One-dimensional leads imply electron-electron interac- 
tion in the system, which is described by Luttinger liq- 
uid model with the dimensionless correlation parameter 
g (0 < g < 1). In the case of strong electron-electron in- 
teraction (g < 1/2), perturbation theory calculations in 



"bare" level width of the fermionic level of QD are valid 
even at low temperatures (i.e. the sequential tunneling 
limit is maintained at all temperatures). This results in 
the strong suppression of tunneling probabilities, due to 
the well-known Kane-Fisher effect^ and in the vanish- 
ing of the conductance at T = 0. This is not the case 
for weak electron-electron interaction (1 > g > 1/2), 
where one could observe a resonant regime of electron 
transport at low temperatures in the case of symmet- 
ric tunnel barrier o 8 ' 10 . The consequence of this fact is 
the perfect (Go = e 2 /h) zero-temperature conductance 
at V -> O 8 -^. 

Thus, it is interesting to understand, how the interplay 
of the effects of strong electron-electron and electron- 
vibron interactions influences electron transport through 
a single-level QD. Although some results in this direc- 
tion have been already obtained for the sequential tun- 
neling limi t 11 ' 12 , these considerations are unable to give 
us correct predictions for the special g — 1/2-case. It was 
shown 8 '! 3 that in the absence of electron- vibron interac- 
tion, a possibility of resonant transport of strongly inter- 
acting electrons (with conductance quantum at T — > 
for the symmetric case) still exists for g > 1/2. At the 
special value g = 1/2 of LL correlation parameter the 
problem of electron transport in such a strongly interact- 
ing system is reduced to the scattering problem for the 
noninteracting fermions, which is exactly solvabl o 13 ' 14 . 
Furthermore, it is impossible to obtain this solution by 
any perturbative method. The origin of this fact lies in 
the special type of symmetry in the correlations between 
the electrons from different leads, which is specific only 
for g = 1/2-system. The generalization of the g = 1/2- 
mode l 13 ' 14 to the case of weak electron- vibron coupling is 
proposed in Ref . [l5[ • The solution of Ref. [ll| treats weak 
electron-vibron coupling pcrturbatively. Thus, it does 
not describe the effects of resonant polaron-assisted tun- 
neling, which take place at sufficiently strong electron- 
vibron coupling. Therefore, a problem of resonant tun- 
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neling for g — 1/2-model with strong electron-vibron in- 
teraction still needs further considerations. 

Below, we consider the influence of quantum vibrations 
on resonant electron transport through a single-level QD, 
weakly coupled to the g = 1/2- Luttinger liquid leads. 
In the polaronic approximation analytical formulae are 
derived for the effective transmission coefficient and dif- 
ferential conductance. The solution describes both the 
resonant tunneling regime for strong electron-vibron in- 
teraction and the sequential tunneling. In particular, we 
have reproduced a known nonperturbative result for elas- 
tic tunneling G(V -> 0, T -> 0) = G = e 2 /h for sym- 
metric junction. Interestingly, this is not the case for 
general (asymmetric) QD, where a Kane-Fisher suppres- 
sion of conductance peak at zero bias still takes place. 
In the strong coupling regime, a correct description of 
all "polaron-assisted" resonances, appearing at energies 
£; = huj l, (I = ±1,2,..) is obtained. At the zero- 
temperature limit, the amplitude of each polron-assisted 
resonant peak in symmetric junction reaches conductance 
quantum, signalling a perfect transmission of polaron at 
resonance energies. We show that the widths of res- 
onances behave non-monotonically, as the functions of 
electron-vibron interaction constant. Moreover, the de- 
tails of such behaviour in our model differ from the pre- 
dictions of a similar model with noninteracting leads 1 
(g = 1 case). Particularly, an additional narrowing 
of widths of all polaron-assisted resonances induced by 
electron-electron interaction is found. A physical in- 
terpretation of such narrowing consists in the special 
type of strong correlations between electrons from dif- 
ferent leads, which affect the tunneling by the additional 
(polaron-assisted) channels. Thus, one can distinguish 
the resonant polaron-assisted tunneling in the interact- 
ing (g = 1/2) and noninteracting (g = 1) systems. 

II. THE MODEL 

The Hamiltonian of our model reads 

H= £ H^+H d + £ 6®. (1) 

j=L,R j=L,R 

In Eq.(l) the first term describes the Hamiltonian of 
one-dimensional infinite Luttinger liquid leads, which is 
quadratic in the bosonized form (we consider the spin- 
less electrons) Hf' = 1/4tt / dx(d x $j(x)) 2 . (Here and 
below we put Hv g = Hvp/g = 1 with "bare" Fermi ve- 
locity vjr.) The bosonic phase fields of the j-th 
lead (j = L,R) are connected with corresponding field 
operators ^j(x) of chiral fermions on the j-th lead in 
strongly nonlinear way by the standard bosonization for- 
mula ^j(x) = exp(i<F,- [x)/sjg) j \f2-Kao (where ag is the 
corresponding lattice constant) 13 . The chiral charge den- 
sity operator pj(x) at the point x of each j-th lead can 
be performed in both fcrmionic and bosonic representa- 
tions as follows p~j(x) = $ff(x)$j(x) — d x ^j(x)/2n^/g. 



Here g - is the dimensionless Luttinger liquid correlation 
parameter. Fermionic field operator $j(x) describes the 
annihilation of chiral fermion, living in one-dimensional 
infinite j-th channel (j = L, R), corresponding to the j-th 
physical lead of the system. However, a negative half-axis 
x < of each j-th channel corresponds to the incoming 
particles moving from the infinity towards the boundary 
(i.e. to the point x = 0~), while the positive one stands 
for the particles, which move back from the boundary (i.e. 
from the point x = + ) to the infinity in the j-th chan- 
nel. The following relation is defined at the boundary 
(in the vicinity of the point x = of each j channel) 13 
#,•(0) = (%(0-) + *j(0+))/2. It is evident from the 
above and due to the continuity of chiral electric charge 
distribution along the system, that the average current 
through the system (i.e. from the left electrode to the 
right one or vice versa) can be performed in a very trans- 
parent form in terms of chiral density operators pj (x) in 
the arbitrary (left or right) leacU^ 

I=(Il(r)) = (e/H){(p L (R)(0-)) - (pl(r)(0 + ))}- (2) 

The second term of Eq.(l) describes the Hamiltonian of 
single-level quantum dot (QD) coupled to a quantum har- 
monic oscillator which models the vibrations of QD. Fol- 
lowing Refs. |l3l l l4j we also take into account Coulomb 
interaction between the leads and QD. 

H d = {Ad+d + Xcd+d *j + ( )*j(°)} + 

j=L,R 

+hu {^(f + x 2 ) + Xxd + d}. (3) 

In Eq.(3) d + (d) - are the standard fermionic creation (an- 
nihilation) operators ({d,d + } = 1) for the electron on 
QD, A - is the level energy, Huiq is the energy of vibra- 
tional quantum (here uiq is the self frequency of quantum 
oscillator) and p,x are the dimensionless bosonic opera- 
tors of the momentum and center-of-mass coordinate of 
QD. The interaction in the QD is controlled by two inde- 
pendent parameters: (i) A is the dimensionless electron- 
vibron coupling constant, and (ii) Ac is the dimension- 
less constant of Coulomb interaction. In our model we 
have fixed Ac at so called Toulouse point (Ac = 2tt), 
while A remains as a free parameter (below we will be 
interested mostly in the case of strong electron-vibron 
coupling A > 1). The last term in Eq.(l) represents the 
tunneling Hamiltonian 

={ lj d+^M + h.c). (4) 

Here jj are the tunneling amplitudes. 

The Hamiltonian H of Eqs.(l-4) can be trans- 
formed into a more convenient form H by intro- 
ducing new symmetric and antisymmetric phase 
fields $±(x) = ($l(x) ± $h(x))/2 and by ap- 
plying of two commuting unitary transforma- 
tions U f = exp[-i{d+d - l/2)$+(0)/ v / 2g] and 
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U b = ex.p(-i\pd+d), H = Uf^HU^JJj 1 . Here 

the transformation removes the electron- vibron 
interaction term from the Hamiltonian (3) of the QD. 
As it was shown 13 , by applying unitary transformation 
Uf to the Hamiltonian (1) at g = 1/2, one can rewrite it 
in terms of new fermions i&±(x) — exp(i$±(x))/ \/2nao 
and remove the <5>+(x)- phase field from the tunneling 
termi 3 -. At the Toulouse point Ac = 2tt one can 
remove also the Coulomb interaction term from the 
transformed Hamiltonian of QD. Thus, if g = 1/2 
and Ac = 27r, one could rewrite the total transformed 

Hamiltonian in the form H = Hi + Hd + H t , where 
Hi = E± l/2vr / dx{d x $±{x)) 2 and 

H d = hd+d+^(f+x*) (5) 

are quadratic now (here A = A — (A 2 /2)Hujq). The trans- 
formed tunneling Hamiltonian takes the form 

H t = d+A+[ 7i *_(0) + 7**1(0)] + 

+ [7L^(0)+7i?*-(0)]Ad. (6) 

Here the operator X = exp(i\p) describes the influ- 
ence of electron-vibron interaction on tunneling, while 
the coupling terms d + X + j R $!±(0) and 7 fl *_(0)A > d 
reveal the existence of additional Andreev-like tun- 
neling of *_-fermions. Operators *±(x) stand for 
new fermions and fulfill standard fermionic anticom- 
mutation relations in the Schroedinger representation 
{^±{x), *±(x )} = S(x — x ). It is evident from the 
above, that new fermionic field operators *±(x)(*jjl(x )) 
annihilate (create) a new nonlocal fermion, which exists 
simultaneously on both (left and right) physical leads 
of the system on the distance |x| from the QD (i.e. 
from the point x = 0) 13 ' 14 . Thus, special electron- 
electron interaction in the system entangles real electrons 
from different (left and right) physical leads in a suffi- 
ciently nonlinear way, making them strongly correlated. 
Density operators p+(x) and p—(x) for i5f± - fermions 
p±(x) = *;£(x)* ± (x) = d x $±{x)/V2~Tr define chiral 
charge- (/>+(x)) and current- (p_ (x)) densities measured 
in the symmetric points on the leads on the distance |x| 
from the QD, p±(x) — Pl(x) ± p R (x). Since $+(x) - 
channel is decoupled now, one can represent the average 
current as follow s 13 ' 14 

(J) = (e/?i){<*±*-(0+)) - (*±*_(0"))}. (7) 



III. QEM-METHOD AND FERMION-BOSON 
FACTORIZATION 

To solve the described model we use a well-known 
quantum equation of motion (QEM) method. The 
Heisenberg equations for the fermionic operators take the 



form 

ihd t d = Ad + l+[ 7i *_(0) + 7**1(0)] (8) 

where, following Refs. |l3l l l4| we defined 
*_(0) = (*-(0") + *_(0+))/2, and 

ihd t ti-(x) = -i9 a *_(x) + 5{x)[ lL Xd- lR X+d+], (9) 

here S(x) is the delta function. Integrating Eq.(9) in the 
vicinity of point x — 0, one obtains 

i[*_(0+) - *_(0")] - -f L Xd- lR X+d+. (10) 

In the absence of electron-vibron coupling 
(X+ =1 = 1), Eqs.(8-10) are reduced to Eq.(6) 
from Ref . [l3| . The formal solution of Eq.(8) can be 
written as follow s 18 ' 19 

d(t) = -(i/h) lim / dt'A + (t')[ 7L #_(0;t') + 

a-tf)J 

+ lR i>t{Q-,t)]e- l{K ~ la){t - t ' )/h (11) 

(with a being positive infinitesimal). Now, substituting 
Eq.(ll) into Eq.(10), we are getting the following basic 
equation for the averages 

fi{($t((r;t)*-(0+;t)) - ($±(0-;f)$_(0-;t))} = 

-Km f <ft'{[7^*l(0-;t)X+(t')i-(t)*_(0;t')> 

+jLlB^t(0-;t)X + (t)X(t)^t(0;t'))}e- l(A - la ^ t - t '^ h 
+h R {^t(0-;t)X+(t)X(t')4>.(0;t)) 
+jLlR(^t(0-;t)X + (t)X(t)^t(0;t))}e^ A+ia ^ t - t '^ h } 

(12) 

here symbol (..) stands for the averaging with the total 
transformed Hamiltonian. 

Now, the central Eq.(12) should be complemented by 
the corresponding equation for bosonic operator p. The 
Heisenberg equation for vibronic subsystem reads 

[d 2 t +iul]p = i^L[X+d+( 7L i>-(0) + 7 «*1(0)) 
-i"d(7£*l(0) + 7fl*-(Q))] 

(13) 

and, obviously, it could be rewritten in the form 

ld?+^]p = ^d t (d+d). (w) 

To proceed further we need to use certain approxima- 
tions. The most evident simplification is to put the right- 
hand side of Eq.(14) to be equal to zero. This corresponds 
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to the case, where bosonic subsystem is effectively unaf- 
fected by fermionic one and the Eq.(14) has a free solu- 
tion 



p (t) = ^(6 + e iuot 



b e 



-iujot\ 



(15) 



Here the operators &J(&o) describe the creation (annihila- 
tion) of a free vibron and fulfill standard bosonic commu- 
tation relation [bo, 6g] = 1. One can see from Eqs. (13, 14), 
that the approximation (15) is always valid in the per- 
turbation theory on T = 7^ + 7# (ro <C huj ,T,eV), 
where the lowest energy scale is Tq -the "bare" width of 
the fermionic level of QD in the "Wide-Band Approxima- 
tion" (WB A- limit). Another approach where one can use 
the Hamiltonian of free vibrons is the so-called "polaron 
tunneling approximation" 17 i 22 . It is valid, when the char- 
acteristic lifetime of electron on QD (~ ft/To) is much 
greater than the time of polaron formation (~ l/A 2 wo). 
In this case only the polaronic states "live" on the dot 
and one can use fermion-boson factorization when cal- 
culating the Green function of polaron. This approach 
allows one to consider resonant tunneling in the system 
with strong electron- vibron interaction (A > 1). Thus, 
when evaluating the averages in Eq.(12) we will assume 
that 



(*±(0-;t)X + (jt )!(t)*_(0;O> ~ 
(*l(0-;t)*_(0;i')). (X + (t')X(t))). 



(16) 



Here the symbols 



and 



stand for the aver- 



^d Hi 

aging with quadratic Hamiltonian of the QD (Hd) and 

of the leads (Hi). In the polaronic approach of Eq.(16) 
it is natural to regard vibronic subsystem as thermally 
equilibrated at the temperature T. 



(btbo) &i =n b (P) = (exp(P) 



1) 



(17) 



with P = hjJo/T. Obviously, our approximations (15- 
17) allow one to generalize the scattering approach elab- 
orated in Ref.[13] for resonant tunneling of interacting 
electrons to the case of resonant tunneling of polarons. 
Under the accepted approximations the solution of ba- 
sic Eq.(12) is formulated in terms of noninteracting - 
fermions backscattered from the point x = or transmit- 
ted through it from negative (x < 0) to positive (x > 0) 
half-axis. Following Ref. 13] we use in Eq.(12) the stan- 
dard momentum decomposition for fermionic field oper- 
ator $-(x,t) 



4>_(x;t) 



dk 
— ( 

2tt 



,ik(t— x) 



Ofe, 

bk, 



x < 
x > 



(18) 



(we set here h = 1 and v g = vp/g = 1) In Eq.(18) a k (dk) 
are the standard fermionic creation (annihilation) oper- 
ators. It is evident, that the average occupation number 
(d k dk) for a new fermionic state at the energy k will be 



as follows (a^a k ')^ = [np(k—fiL)—np(k—fiR)]S(k—k ), 

where n F (k-^ L{R) ) = (exp(/3(/s - Mi(fl))) + are the 
Fermi distribution functions of electron in the L(R) reser- 
voirs. Analogously to Ref.[l3], we represent operator bk 
in Eq.(18) as bk = t(k)dk, where t(k) is the transmis- 
sion amplitude. The reflection amplitude r(k) is defined 
by the equation a^ k = r(k)dk and describes the process 
of the Andreev-like reflection of the incoming quasiparti- 
cle. We put also b_ k = 0, since there are no transmitted 
"holes" in the system. The considerations given above 
allow us to rewrite a basic formula (7) for the average 
current through our system. It takes the Landauer-type 
formic 

I(V) = (I) = ^J deR(e)[n F (e - eV) - n F (e)], (19) 

where R(s) = 1 — |t(e)| 2 is the energy-dependent re- 
flection coefficient for fermions, which determines 
the effective transmission coefficient for physical electrons 
transferred through the QD. 



IV. RESULTS AND DISCUSSION 

One can see that the problem of electron transport in 
the considered model is reduced now to the evaluation of 
the effective transmission coefficient R(e). Thus, regard- 
ing the basic Eq.(12) with its complex conjugated equa- 
tion, after the averaging under the assumption about a 
fermion-boson factorization (15)-(18), one can solve the 
resulting integral equation with respect to the complex 
transmission amplitude t(e) of fermions. Partic- 
ularly, after the integrating over dt , taking the limits 
t — >• oo and a — >• 0, it is possible to derive the follow- 
ing basic formula for the effective transmission coefficient 
R(e) from Eq.(19) 



R(s) = 



4r 2 



[(1 + B\ )(1 + B 2 _) + 2r 2 (i + B+B_) + r 4 ] 



(20) 



Here 



B± = B ± (e,P)= 



T Fi(P) ll± (s) 



A? 



+ 00 



T = f(e,p)= 



f=-c 



r F t (^) 7 (e) 
A 2 -e 2 



(21) 



with 



7/±(e) 



A,±e 
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lie) 



1L1R . . 

.,2 i -,2 I fc ' 



7£ + 7fl 



(22) 



and 



F,C8) - e- A2(1+2 " 6) / i (AV^(l + ^))e- W2 . (23) 

In Eqs.(21,23) Ii{z) is the Bessel function of Z-th order of 
the imaginary argument. In Eqs.(22) A/ = A — fnoql is 
the energy of the Z-th resonance. 

Formulas (20)-(23) generalize the results obtained by 
Komnik and Gogolin in Ref.[l3T| on the case of strong 
electron-vibron interaction in the QD and represent the 
basic result of this paper. Indeed, in the case when 
A = (i.e. in the absence of electron-vibron interaction) 
our general formula (20) is reduced to the basic result 
(Eq.(10)) of Ref.[l3| for the general case of asymmetric 
tunnel barriers (71, ^ Jr). One can see that the influ- 
ence of the asymmetry of tunnel barriers on the trans- 
port properties of the system is concerned mostly in the 
renormalization ~ r/Rs{e) of the effective transmission 
coefficient i?s(e) for the symmetric junction (7^ = 7^). 
Here 77 < 1 is the asymmetry parameter 77 — (jl/ir) 2 if 

7l < 7R) and V = {ir/il) 2 if 7R < 1l- 

In the asymmetric case, as one can see from Fig. la, the 
"zero-bias" (with 1 = 0) resonance has a dip at e — s- 0, 
if A = 0. This dip shrinks to zero, when e — at the 
arbitrary value of asymmetry parameter 77 < 1. The 
presence of such dip at e — > in the case of even small 
asymmetry of tunnel barriers is the manifestation of the 
Kane-Fisher effecl£i2 at e — > ef = 0. This is because at 
77 ^ 1 the specific particle-hole symmetry in our model 
breaks down, and more common Luttinger liquid physics 
is revealed. 

Thus, all distinguishing features of the model can be 
obtained by considering a more simple case of symmetric 
tunnel barriers 7^ = jr. In this case our central expres- 
sion (20) for the effective transmission coefficient i?s(e) 
takes a very transparent form 



Rs(e) = 



1 



1 + (S(e)Y 



with 



S(e) 



g £ r F ; (/3) 

l= — oc 



A? 



(24) 



(25) 



Evidently, in the absence of electron-vibron interac- 
tion (A = 0), formulae (24,25) reproduce the result for 
resonant tunneling between Luttinger liquid leads with 
g = 1/2 in the symmetric cas o 10 i 13 ' 14 . Particularly, at 
A = and A = 0, Eq.(24) turns into the usual Breit- 
Wigner expression for transmission coefficient of nonin- 
teracting electrons. This is because, in that case the 
mapping on the model with Fermi leads (g — 1) becomes 
valio~. 

Now let us analyze formulae (24), (25) in a more general 
situation, where A 7^ 0. It is worth to point out, that by 



R(e) 



a) 



R(e) 




FIG. 1: The effective transmission coefficient R(e) of Eq.(20) 
as the function of energy (in the units of ftuo), for differ- 
ent values of the asymmetry parameter r\ = T h ( R )/T R ( L }, in 
comparison with effective transmission coefficient (red dotted 
line) obtained in Ref.[H (Eqs. (22,23)) for the g = 1- model 
with symmetric tunnel barriers (77 = 1). On Fig. la zero-bias 
peak is depicted for: 77 = 1 (symmetric case) - blue solid line; 
r\ = 0.4 - blue dotted line; r\ = 0.1 - magenta dotted line; 
and r? = 0.025 (strongly asymmetric case) - black solid line. 
On Fig. lb, in addition to the zero- bias peak, two "polaron- 
assisted" peaks are shown. These peaks describe resonant 
Andreev-like polaron-assisted tunneling. Here r\ = 1 (sym- 
metric case) - blue solid line; 77 = 0.1 - magenta dotted line; 
and rj = 0.025 (strongly asymmetric case) - black solid line. 
On Figs.la,b we also put A 2 = 1, To/fii^o = 0.25. 



means of the physically transparent method formulated 
above, one could also solve the model with Fermi- leads 
(g = 1- case) for strong electron-vibron interaction A > 1 
(see Ref.[l7j])- In particular, for the symmetric junction 
{lL = Ir = 7o)j eliminating the exotic coupling terms 
(l + X + 7o*l (0) and 70*- (0) X d from the tunnel Hamilto- 
nian (6), and applying the entire scheme described above, 
one could reproduce the central results (Eqs. (22), (23)) 
of Ref.[17], obtained there by means of more rigorous 
mathematical methods (Full Counting Statistics (FCS) 
and Keldysh technique). This is because, the central 
approach of Ref.[l7j (the linear diagram resummation 
scheme) is included into our polaronic approximation of 
Eqs.(15-18). 

One can see, that in general case (if A 7^ and A/ 7^ 0), 
the expressions (24-25) differ from corresponding formu- 
lae (Eqs. (22-23)) of Ref.0] for the case of Fermi leads 
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(with g — 1). Indeed, in our model a special symme- 
try: Pl(b,)( x ) ~Pr(l){ x ) °f charge density excita- 
tions on different physical leads takes place. Such sym- 
metry makes electrons from different leads be strongly 
correlated in the vicinity of QD and implies a special 
(Andreev-like) type of resonant tunneling in the system 
at A; ^ and jl = 1r- In the latter perfect 
transmission of the electron through QD, along the in- 
finite one-dimensional system (when in Eq.(19) one has 
i?s(e) — 1)) can be treated as the total Andreev reflec- 
tion of spatially nonlocal fermion from the bound- 
ary at the point x = back to the half-infinite ID 
system with x < 0. During such a reflection the in- 
coming particle-type excitations are transformed into the 
opposite-moving excitations of the hole-type at the same 
energies. 

The energy- dependent effective transmission coeffi- 
cient of Eq.(20) is plotted on Fig.l for different values 
of the asymmetry parameter r\ = T T i n\ /^r(l) i m com- 
parison with one, calculated in Ref.[l7| for the case of 
Fermi leads {g — 1) and symmetric tunnel barriers. It is 
essential, that in the resonant tunneling limit at T = all 
obtained from Eqs. (20,24) resonant peaks describe elas- 
tic processes of a perfect transmission of the polaron at 
resonant energies ei = HujqI, I = 0, ±1,2, ... In such 
processes, one or more virtual vibrons are emitted and, 
then absorbed, leaving the fermionic subsystem of the 
QD in the same quantum state as before the polaron 
transmission 17 . On the other hand, in the opposite case 
of sequential tunneling (Tlh <C T <C Huio) only the in- 
elastic tunneling, accompanied by the emission of real 
vibrons, produce small (satellite) resonances, while all 
virtual processes result only in the "polaronic narrow- 
ing" of all resonances^. 

Regarding the off-resonant energies: 

hu>o ^> \si — A; | 3> To, one may conclude from 
Fig. lb, that the transmission coefficient (20) shrinks 
to zero at these energies if A 2 To 3> T — > 0. This is 
the consequence of destructive interference of different 
virtual polaronic states on the QDi 7 -. As one can see 
from Fig.l, although a perfect transmission at zero 
temperature in both cases (for g = 1/2- and g = 1- 
models) takes place at the same energies ei = HuiqI, 
I = 0, ±1,2, .., the electron-electron interaction in the 
g = 1/2- case sufficiently narrows all "polaron-assisted" 
(with I > 1) resonances, as compared to the case of Fermi 
leads (g = 1). This fact clearly shows, that although a 
resonant polaron-assisted tunneling is possible for both 
systems with g = 1- and g = 1/2- leads, the details of 
such processes are different in these two cases. Indeed, 
the resonant tunneling in g = 1/2- case requires from 
the electrons placed on different leads to be in definite 
strongly correlated quantum states during the process of 
resonant tunneling. These strong correlations entangle 
physical electrons from different leads in the vicinity of 
QD. Obviously, the probability of resonant state, which 
involves more than one physical electron (when g = 1/2) 
at resonant energy far from ep = (i.e. in the case I > 1) 



is expected to be much smaller than the probability of 
the resonant state, which involves only a single electron 
of the same energy (when g = 1). To estimate this effect 
quantitatively for the case Tq -C Hloq -C eV, one could 
perform the effective transmission coefficient of Eq.(24) 
at resonant energy ei near the l-th polaron-assisted 
resonance (/ > 1) as follows 

Rso(ei) * =5 T l f !l ] -y 2 (26) 

T eff(l) + ( A * " £ <) 

with the following approximate expression for the effec- 
tive width T e ffm of each Z-th (I = 1, 2, .., l m ) resonance 

Tefm ^ i + t'aW (27) 

Here. 

r ; = r oe - A2 (A 2 )7/! (28) 

-is a low-temperature asymptotic for the effective width 
of the l-th resonance in the case of Fermi leads^ 7 - and 
A/ = HujqI, (A = 0), < I < l m , with the maximal num- 
ber of vibrons emitted l m ~ [eV/huj ). This is because, 
at T, r <C there are no processes of absorption of 
any "external" vibrons and since the energies of incoming 
quasiparticles are limited by eV at T <C eV. 

The important consequence of Eqs. (27,28) is that, the 
obtained relative additional narrowing T e ff^/Ti due 
to electron-electron interaction is found to be strong 
(~ 1/2) and approximately the same for all polaron- 
assisted resonances. This feature seems to be specific 
only for the considered model and might help us in the 
distinguishing of this important g = 1/2- case among 
other experimental realizations of ID molecular transis- 
tors. This is clearly illustrated by Fig. 2, where the effec- 
tive widths of polaron-assisted resonances are plotted (in 
the units of To), as functions of electron- vibron interac- 
tion A 2 for different values of I > 1 in both g = 1 /2- and 
g = 1- cases. 

From Fig. 2 one can easily observe two main effects. 
First of all, one can see the described above effect of the 
additional narrowing of all polaron-assisted resonances, 
due to specific electron-electron interaction in the sys- 
tem with g = 1/2- LL leads. As it could be seen from 
Eq.(27), this effect is strong enough at I > 1. Besides 
that, for strong electron- vibron interaction (if A > 1) 
in both g = 1/2- and g = 1- cases a strongly non- 
monotonic behaviour of the effective resonance widths 
(with the maximum at A 2 ~ I) as functions of A 2 takes 
place. The latter phenomenon is a well-known conse- 
quence of the fermion-boson factorization procedure and 
is caused by the interplay of the effects of vibron-assisted 
tunneling and polaronic blockad o 1 ^ 1 . Although our gen- 
eral formulae are valid in the limiting case A = for both 
resonant and off-resonant tunneling, for A — > (A / 0), 
the predictions of our method coincide with the results 
of perturbation theory calculations in small A 2 for each 
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FIG. 2: The effective widths of different polaron-assisted res- 
onances from Eq.(27)(in the units of To) for symmetric junc- 
tion (jj = 1), as the functions of electron- vibron interaction 
A 2 . Blue solid line corresponds to effective width of polaron- 
assisted resonance for g = 1/2-model with 1 = 1, while the 
brown solid line represents the resonance with I — 3 for the 
same model. Red dotted line and magenta dotted line cor- 
respond to widths of polaron-assisted resonances with I = 1 
and I = 3, obtained for the g = 1-modei^. (Here we put 

r /fiw = o.i.) 



I ^ only while T/T > (\ 2 /l) 1 . Thus, it is reason- 
able to conclude, that if A <C 1 (but A ^ 0) the case 
T/T < (\ 2 /l) 1 < 1 (A -> 0, T = 0) is not described by 
means of our method^. 

Finally, let us consider the behaviour of the differential 
conductance Gso(V) = dI(V)/dV for our model in the 
simplest symmetric case (jl = 7_r)- At high tempera- 
tures X 2 hu>Q <C T, when all thermally activated (vibronic) 
channels contribute to electron tunneling, the polaronic 
blockade in the zero-bias peak is totally lifte d 20 ' 21 and we 
get a standard high-temperature asymptotic ~ GqTq/T 
for the conductance of the g — 1-modcl with noninteract- 
ing electrons (here Go = e 2 /h is the conductance quan- 
tum). This is because, at high temperatures all specific 
quantum features of both electron-electron and electron- 
vibron interactions are " smeared out" by thermal fluctu- 
ations. Much more interesting is the case of low tempera- 
tures, where T <C hujQ. In this case, a main contribution 
to effective transmission coefficient (24) goes from the 
resonant terms of the form (26). Therefore, one could 
obtain a following low-temperature asymptotic for differ- 
ential conductance Gso(V) 



level (I = 0, 1,2, ..,l m ) (-here the "detuning" term A is 
insufficient, due to the existence of the gate voltage, thus 
we put it to be equal to zero). The effective widths ^ e ff(i) 
are from Eq.(27), and the maximal number l m of vibrons 
emitted is of the order of the integer part of quantity 
[eV/Hujo]. The resulting differential conductance (29), 
as the function of bias V, exhibits a sequence of sharp 
resonances at eV = eVt m — A; m (l m — 0,1,2,..). The 
main contribution to the sum in Eq.(29) goes from the 
resonant term with l m ~ [ey/fiwo] if T -C hcu <C eV, 
and from the term with l m = if T -C eV -C fv^o, 
correspondingly. Thus, differential conductance (29), as 
the function of resonant values of bias eVi m — A/ m can 
be estimated for these two cases as follows 



tanh 



4T C 



-| In 



Gso(Vi m ) « G < 



tanh 



r e" 



17 l m =0 
^ (T < eVi m < hu ) 



Evidently, 
(where T 



(30) 

the Eq.(30) describes both resonant 



eff(lv 



and sequential (where 



T F e ffn\) regimes of polaron tunneling as well. 
Indeed, if eVi m = eVb <C fu^o (l m — 0) then 
Gso(Vb) ~ Gor e~ A "/T. This formula reproduces a 
well-know n 6 ' 7 ' 11 ' 16 r^/T-scaling for the low-temperature 
conductance of the g = 1-model, with the renormalized 
bare level width T\ — Toe" , as the result of polaronic 
(Frank-Condon) blockade in the elastic channel of tun- 
neling. For differential conductance (eVi m 3> ftwo) i n the 
sequential tunneling limit T > (l/2)(A 2 /Z m ) /m r e -A2 , 
(l m ~ [eVim/hoJo] ^> 1 and T <C fiwo <C eVJ m ) the situa- 
tion is quite similar. The expected 1/T-scaling of differ- 
ential conductance (30) in this case will be the following 



G SO (Vim) 



— A 2 



T 



X 2 Hujq 

eVi m 



eVim — huj()l ri 



(31) 



G5o(^) = G ^i{tanh 



1=0 



tanh 



Formula (31) gives us the "hights" of low-temperature 
differential conductance peaks in the sequential tunnel- 
ing limit, at resonant values of bias: eVi m ~ fiwo'm, 
l m = 1,2,.. in the case of symmetric junction. The 
only difference between expression (31) and correspond- 
ing limiting case for the g = 1- model is a pre-factor 
~ 1/2 in Eq.(31). It results in a sufficient suppression 
of all satellite (with l m > 1) peaks, as compared to the 
In Eq.(29) A/ = HujqI is the energy of the l-th resonant g — 1- case. Naturally, the above expansion of Eq.(30) 



A t -eV + r e//(0 /2 
2T 

&i-eV-T efm /2 
2T 



(29) 
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eV 

hujQ 

FIG. 3: The low-temperature differential conductance 
Gso(V) (in the units of conductance quantum Go) from 
Eq.(29), as the function of bias V (in the units of hwo/e), 
for different values of electron- vibron interaction A 2 . Here the 
blue solid line represents the case A 2 = 4, while the red dashed 
line corresponds to A 2 = 6, and the brown solid line describes 
the case A 2 = 8 (-strong electron- vibron interaction) . Also we 
put here friuo/T = 8 and Yo/f\Uo = 0.1. 



for the sequential tunneling limit coincides (up to a pre- 
factor ~ 1/2) with the result of perturbation theory for 
the t/ = l- model in small To (Tq <C T), or in small A 2 
at non-zero temperatures (if A 2 <C T/r )££ Naturally, 
this fact is because of the mapping on the g = 1- model. 
A novel feature here is that, such mapping is not com- 
plete any more for energies e ~ eV > ftwcb even in the 
sequential tunneling limit. Only for eV <C Swo the cor- 
respondence between the g = 1/2- and g = 1- models 
with strong electron-vibron interaction is complete. Of 
course, the standard perturbation theory in small Tq for 
the SET weakly coupled to the LL-leads with arbitrary 
g ?i n i 12 (0 < g < 1) is unable to give a correct description 
for our case, since it does not take into account strong 
correlations between tunneling events, which are specific 
for g = 1/2- system. 

As regards to the resonant tunneling limit 



T < r e//(;m) in Eqs. (29-30), the interplay of the 
effects of electron-electron correlations and strong 
electron-vibron interaction becomes nonperturbativc 
in this case. In particular, in this limit there is no 
significant difference between the peaks of differential 
conductance, corresponded to polaron-assisted (with 
lm 7^ 0) and zero-bias (with l m = 0) channels of 
tunneling. Besides that, strong interactions of both 
types (electron-electron and electron-vibron) lead to a 
sufficient narrowing of all polaron-assisted (with I ^ 0) 
resonances at: A 2 > 1, whereas in the zero-bias channel 
(with I = 0) only a polaronic narrowing takes place 
at: A > 1. In the zero-temperature limit T — > 0, one 
will have a perfect resonant tunneling of a polaron with 
conductance quantum Gso(V; m ) ~^ Gq at resonance 
values of bias eV\ m — A; m , (l m = 1,2,..), but only 



if < T/T < (A 2 /Z m )S (T 0). Thus, evidently, 
the zero-temperature perturbation theory in small A 2 
(see Ref.pjl) does not describe the considered case of 
resonant polaron-assisted tunneling in the g = 1/2- 
model at T = 0. 

On the other hand, in the case of strong electron- 
vibron interaction (A > 1) at low, but nonzero tem- 
peratures (0 < T/Tq <C 1) as well, as in the discussed 
above sequential tunneling limit, one can observe in the 
model the domination of vibron-assisted electron trans- 
port (see Eqs. (29-30) and Fig. 3). This effect was pre- 
dicted earlie r 11 ' 12 , but only in the limits of perturbation 
theory in Tq. On Fig. 3 the differential conductance of 
Eq.(29) is plotted (in the units of Go), as the function 
of bias voltage V for the "intermediate" region of tem- 
peratures, between the "resonant" and the "sequential" 
tunneling regimes (when T/Tq > 1), in the case of strong 
electron-vibron interaction (A > 1). One can see from 
Fig. 3, that polaron-assisted (with / ^ 0) resonant peaks 
strongly dominate the zero-bias (with I = 0) peak for all 
values of A > 1. Therefore, the highest peak on Fig. 3 
is always polaron-assisted and corresponds to the case 
where eV = eVi m ~ X 2 hidQ, for every given value of A 2 
(A>1). 

In the zero-temperature limit, the resonant average 
current through the system will be the following 



T (Vi m ) w r r e- A2 

h 



1=0 



l\ 



l m ~[eV lm /Hu> ] 



(32) 

with the bias-dependent effective width of the highest 
virtual resonant level. It is evident from Eq.(32), that 
zero-temperature average current as the function of bias 
reaches its maximal value Iomax = (e/h)To in the limit 
eV/huio — > oo, when all possible bias-activated virtual 
channels of polaron-assisted tunneling are opened and, 
as a result, the polaronic blockade is totally lifte d 20 ' 21 . 
Thus, it is reasonable to conclude that, in the case of 
symmetric junction, the additional narrowing of Eq.(27), 
as well as the existence of perfect transmission of a po- 
laron at sequence of resonance energies hwol, I = 0,1,2,.. 
in the zero-temperature limit (T — > 0), may serve as the 
manifestation of the novel (Andreev-like) type of polaron- 
assisted resonant tunneling of strongly interacting elec- 
trons. We think, that the resonant tunneling of such 
a type is unique for the considered g = 1/2-model with 
strong electron-vibron interaction and represents a conse- 
quence of strong correlations between the electrons from 
different physical leads of the system. 



V. SUMMARY 

In the above, the resonant tunneling of strongly inter- 
acting electrons through a single-level vibrating quantum 
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dot (QD) is considered for the case of strong electron- 
vibron interaction in the QD. Corresponding transport 
problem is solved in terms of the scattering of noninter- 
acting fermions, which involve the "entangled" electrons 
from both physical leads of the system in the vicinity of 
QD. As a result, the general formulae for the effective 
transmission coefficient and for differential conductance 
are obtained. It is found, that in the case of symmetric 
junction, for strong electron-vibron interaction, at suffi- 
ciently low and zero temperatures, a novel (Andreev-like) 
type of resonant polaron-assisted tunneling is realized. It 
turns out, that some features of this type of tunneling are 
quite similar to ones for the noninteracting electrons (in 
the case of Fermi leads, where g = 1). Especially, such 
resonant tunneling is characterized at zero temperature 
by perfect transmission (with conductance quantum) of a 
polaron at sequence of resonance voltages eVi m « hw^lm, 
l m = 0, 1, 2, ... The effective widths of all polaron-assisted 
(with l m ^ 0) resonances depend non-monotonically on 
electron-vibron interaction constant A, with the maxi- 
mum at A 2 ~ l m for each / m -th (l m ^ 0) resonance. 
This feature leads to the domination of polaron-assisted 
electron transport in the case of strong electron-vibron 
interaction. But, despite partial mapping on the non- 



interacting (g = 1) case, the most important difference 
between resonant tunneling in both g — 1/2- and g = 1- 
models is the additional narrowing of the widths of all 
polaron-assisted resonances in the g — 1/2- model, as 
compared to the g = 1- case. This relative additional 
narrowing is found to be strong and roughly the same 
(~ 1/2) for all polaron-assisted resonances with: l m > 1. 

Such novel feature points out on special mechanism 
of polaron-assisted resonant tunneling, which seems to 
be unique for the considered g — 1/2- model. Partic- 
ularly, in the case of g = 1/2- LL leads, the resonant 
quantum state of certain energy represents the "entan- 
gled" quantum state of the electrons from both physical 
leads of the system. As the consequence, physical elec- 
trons from different leads of the system become strongly 
correlated, due to the special type of electron-electron 
interaction. It is relevant, that the revealed additional 
narrowing (roughly in 1/2 times) of all polaron-assisted 
resonances, as compared to the case of Fermi leads, may 
serve, as the important distinguishing feature of such spe- 
cial g = 1/2-type of ID electron systems. 

In conclusion, the author would like to thank 
I.V.Krive, and also to S.I.Kulinich, R.I.Shekhter and 
L.A.Pastur for the valuable discussions. 
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